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Abstract 

The impact of the third (skewness) and fourth (kurtosis) reduced 
centered moments on the statistical modehng of El hues in complex 
atomic spectra is investigated through the use of Gram-Charlier, Nor- 
mal Inverse Gaussian and Generalized Gaussian distributions. It is 
shown that the modeling of unresolved transition arrays with non- 
Gaussian distributions may reveal more detailed structures, due essen- 
tially to the large value of the kurtosis. In the present work, focus is put 
essentially on the Generalized Gaussian, the power of the argument in 
the exponential being constrained by the kurtosis value. The relevance 
of the new statistical line distribution is checked by comparisons with 
smoothed detailed line-by-line calculations and through the analysis 
of 2p 3d transitions of recent laser or Z-pinch absorption measure- 
ments. The issue of calculating high-order moments is also discussed 
(Racah algebra, Jucys graphical method, semi-empirical approach ...). 

1 Introduction 

The detailed calculation of all the electric-dipole (El) line energies and ra- 
diative strengths in complex atomic spectra is an overwhelming task. In 
some circumstances, it can even be useless. Indeed, when the density is 
sufficiently high so that the physical broadening mechanisms (Stark, ...) 
are important and/or when the number of lines in an energy range be- 
comes large, the lines coalesce into broad structures. Statistical methods 
are required because, experimentally, some quantities cannot be determined 
individually, but only as weighted average quantities (this is the case for in- 
stance in emission/absorption spectra of highly ionized atoms). Moreover, 
explicit quantum calculations can be inappropriate, e.g. if the eigenvalues of 
the hamiltonian are not known with sufficient precision. In addition, global 
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methods can reveal physical properties hidden by a detailed treatment of 
levels and lines ("one cannot see the wood for the trees") [1]. 

A transition array [2] of El lines is characterized by a specific distribution 
of photon energy E: 

I{E) = A{E) ^(S), (1) 

where the function A{E) assumes that each line is represented by a Dirac 
(5-function: 

A[E) = C,Y.^ {Eabf Sab SiE - E,b) (2) 

a,b 

which, using the appropriate constant factor Cp, represents either the 
opacity (p = 1) or the emissivity of the source {p = 4). The sum runs over 
the upper and lower levels of each line belonging to the transition array. 
The density of ions excited in level a is noted Na, and Qa is the degeneracy 
of level a. The energy of the line a ^ 6 is 

Eab = Eb-Ea =< b\H\b > - < a\H\a >, (3) 

where H is the Hamiltonian of the system and the line strength Sab is 
equal to | < a\Z\b > p, where Z is the z component of the dipole transi- 
tion operator. The normalized profile ^(E — Eab) takes into account the 
broadening of the line in the plasma due to natural width, Doppler effect, 
ionic Stark effect, electron collisions, etc. In the UTA (Unresolved Transi- 
tion Arrays) modeling [3], the discrete distribution A{E) can be replaced 
by a continuous function (usually Gaussian) which preserves its first- and 
second-order moments. The density Na is assumed to be proportional to 
its statistical weight g^: Na « Nga/g, with N = Y^a^a and g = Y.a9a- 
In order to avoid the sum over the term {EabY in Eq- ([21) j the line energy 
Eab is replaced by the center of gravity Eq of the transition array, i. e. 
{EabY ~ {EgY ■ These assumptions allow one to express the moments of 
this distribution as 



A{E) E- dE Ea,b Sab {EabY 
fZoME)dE ^ Y:a,bSab 

Za,b[< b\H\b >-< a\H\a >]"| < a\Z\b > p 
Ea,b\ < ^\2\b > \' ■ 

It is possible to derive analytical formulae for the moments //n(^) using 
Racah's quantum-mechanical algebra and second-quantization techniques of 
Judd (5- Such expressions, which depend only on radial integrals, have been 
published by Bauche-Arnoult, et al. [3l O [6l [T] for the moments fin (with 



MA) = 



2 



n ^ 3) of several kinds of transition arrays (relativistic or not). It is useful 
to introduce the reduced centered moments of the distribution defined by 

an{A) = -j A{E) [^-^) dE, (5) 

where ni is the center-of-gravity of the strength-weighted line energies, 
a = iJ,2 — t^i is the standard deviation and fl = A{E) dE is the total 
area of the distribution. The use of an{A) instead of /U„(A) allows one to 
avoid numerical problems due to the occurence of large numbers. The first 
values are ao = 1) c^i = and 02 = 1- The distribution A[E) is therefore 
fully characterized by the values of 0, /xi, a and of the high-order moments 
an with n > 2. The first four moments are often sufficient to capture 
the global shape of the distribution A{E). The third- and fourth-order re- 
duced centered moments 03 and 04 are named skewness and kurtosis. They 
quantify respectively the asymmetry and sharpness of the distribution. The 
kurtosis is usually compared to the value 04 = 3 for a Gaussian. However, 
the choice of the distribution does not play any role in the derivation of the 
moments. 



2 Impact of the third- and fourth-order moments 
2.1 Gram-Charlier expansion series 

In order to investigate the impact of skewness and kurtosis, it is possible to 
use the Gram-Charlier expansion series |S]: 

n e--2 ' ^ 



GCniE) = - ^ ( 1 + ^ cfc Hefc(u) ] , (6) 



k=2 



with 



mt(fc/2) . 

where u = {E — ^i)/a^ n is the order of the moment and Ilefc(ti) is the 
Hermite polynomial of order k obeying the recursion relation: 

He„+i(x) = X Hen(x) - n He„_i(x), (8) 

with Heo(x) = 1 and Hei(x) = x. The Gram-Charlier expansion series 
uses the reduced centered moments ak{A) of the discrete distribution A{E). 
When carried on to infinity, it can be shown that the Gram-Charlier series is 
an exact representation of the distribution. Since IIe2(x) = x^ — 1, IIe3(x) = 
x^ — 3x and He4(x) = x^ — 6x^ -|- 3, one finds that the fourth-order Gram- 
Charlier series reads: 
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n e-- 




(«4 - 3) 
24 



(S-e-u^ + n^) . (9) 



One of the main disadvantages of the Gram-Charher profile is that in 
certain circumstances (see Fig. [T]), it exhibits some negative features. 

The levels of an electronic configuration verify approximate symmetries, 
which are called couplings in atomic spectroscopy. In fact, highly asymmet- 
rical line distributions can be found due to some large exchange Slater 
integral. This may happen in arrays of the type — l^l', which are very 
numerous, because in highly charged ions, the transitions decaying to the 
ground configuration belong to arrays of this type. The skewness increases 
steeply as a function of atomic number Z along an iso-electronic sequence. 
Assuming hydrogenic behaviour, Slater and spin-orbit integrals vary as Z 
and Z^, respectively. It was shown in that case by Bauche et al. [9j that 03 
behaves roughly as Z^. Therefore, asymmetry is particularly pronounced 
for highly ionized heavy atoms. They can also occur in situations where 
the spin-orbit interactions are strong enough to split the transition array 
(medium- or high-Z elements). In the latter case, the asymmetrical shape 
of the array can also be restored by considering the superposition of sym- 
metrical subarrays. In fact, the impact of skewness is usually small in the 
conditions typical of the laser experiments conducted to date (density of the 
order of 0.01 g/cm^ and temperature of a few tens of eV). Fig. [2] displays 
the comparison between a skewed Gaussian (fourth-order Gram-Charlier 
(see Eq. ([9])) with = 0.6 and 04 = 3 compared to a Gaussian modeling of 
lines in the case of a bromine plasma at T=47 eV and /9=0.04 g/cm^. The 
spectra are almost indistinguishable. 

2.2 Normal Inverse Gaussian 

If one wants to take into account the effect of asymmetry, it can be fruitfully 
done using the NIG (Normal Inverse Gaussian) distribution [2S)l- 



where Ki is a modified Bessel function of the third kind. The four 
parameters a, (3, 5 and fi are obtained directly from the knowledge of fii, 
a, 03 and 04. The corresponding relations as well as the role of the latter 
parameters can be found in Table [TJ Unlike the Gram-Charlier expansion 
series, the NIG (see Fig. [3]) cannot have negative values. Moreover, the 
validity domain of the NIG is wider than the positivity domain of Gram- 
Charlier distribution, as shown in Fig. HI The NIG distribution enables one 
to model symmetric or asymmetric distributions with possibly long tails in 




(10) 
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both directions. The tails are much more prominent than in the Gaussian 
distribution. 




2.3 Generalized Gaussian 

An interesting choice to study solely the effects of the kurtosis is the gener- 
alized Gaussian distribution (GG), defined by: 

GG{E) = —-^^—- with A = W^, (11) 
a 2X r(i) ^' ^'^^ 

where is a positive real number, and T{x) is the ordinary gamma 
function. The even-order moments of a GG function read: 

«2.(GG) = (12) 

whereas the odd-order moments are null, a2k+i{P) = 0, since the GG is 
symmetric. The parameter can be obtained by constraining the kurtosis 
coefficient, and thus solving the equation 

r(i/z.)r(5/z.) 
= r(3/.)2 • 

This distribution is represented in Fig. [S] for several values of the pa- 
rameter 1/ (i.e. for different values of the kurtosis a^). The GG function 
has interesting properties. It is a simple increasing (decreasing) function for 
u < {u > 0), without negative values in contrast with GG series. The 
Gaussian (i/ = 2) and the Laplace {ly = 1) distributions are special cases of 
GG functions with a kurtosis coefficient equal to 3 and 6, respectively. The 
discontinuity of the derivative at u = for < ^ 1 disappears with the 
convolution by another function. The full width at half maximum (FWHM) 
of a GG function is 2 aX (log 2)^ , and therefore depends on the distribution 
itself. For example, the above formula gives FWHM = 2.35 a for a Gaussian 
(z/ = 2) and FWHM = 0.98 a for a Laplace distribution (z^ = 1). The root 
of Eq. ()13p is fairly well approximated by the fitting function 



0.796349 
2 



z/ = 1.62796 [ln(a4 - 0.783143)] "i . (14) 



2.4 Comparisons of the distributions 

Fig. [6] displays the comparison of a transmission spectrum (for the same 
case as Fig. [2]) with GG and GC profiles with a3=0 and 04=6. We observe 
that both distributions give rise to slightly different spectra, although they 
are characterized by the same first four moments. The differences are the 
signature of high-order moments a„ (n > 4). Indeed, the first constained 
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moments fix ipso facto tlie liigher-order unconstrained moments. Tliis can 
be seen in Table [2l for the transition array 3(P — 3d^4p in Br XII, which 
shows the values of Ok up to the order 14 for several distributions. The 
exact values of the moments are calculated with Cowan's atomic structure 
code (routines RCN/RCN2/RCG). They are compared with the values 
for the GG function, using Eq. (I12|) . with the values for the fourth-order 
GC series (obtained by setting = for k > 6) and for the NIG function 
(requiring a = 0.82, (3 = 0, 6 = 3.17 and fj, = 154.18). Functions GG, 
GC4 and NIG are constrained by the moments of order A; < 4 of the exact 
distribution. Values for the Gaussian function are also shown. It is observed, 
for this particular case, that the GG function is slightly closer to the exact 
values than the other distributions. 

Transition arrays of the kind l^l' — l^l" are known to be much sharper 
than a Gaussian, because of the strong selection rules on the core Fig. 
[7] displays a comparison of the line distribution of transition array 
3d^4d for V II calculated with Cowan's code, and modeled by Gaussian, GG, 
GC4 and NIG profiles. All the distributions are convolved by a Gaussian 
having a FWIIM=0.2 eV. We can see that the width and heigth of the 
distribution are not well depicted by the Gaussian. The Gram-Charlier 
GC4 distribution has negative features. The GG and NIG distributions are 
more suited in that case. Fig. [8] displays the same kind of comparisons for 
the transition array 3d^4s ?>d!^Ap of Co V with a Gaussian of FWHM=0.1 
eV. The Gram-Charlier GC4 has not been represented since it exhibits huge 
non-physical bumps. In that case also the GG and NIG distributions give a 
better agreement with the exact distribution than the Gaussian. This visual 
agreement has been confirmed by RMS calculations. Finally, Fig. [9] displays 
a comparison of the profiles for the transition array 3(i'^4p — > "id? Ad for Co VI, 
the distributions being convolved by a Gaussian having a FWHM=0.3 eV. 
In that case the distribution is asymmetric, and therefore the NIG appears 
to be the best distribution, but the GG still gives better results that the 
Gaussian. 

We may conclude that the Gaussian is not a good representation for 
most transition arrays and that GG or NIG are usually better choices than 
GC when accounting for skewness and/or kurtosis coefficients. 

3 Interpretation of recent 2p 3d absorption ex- 
periments 

Chenais-Popovics et al. [12j measured the absorption of the 2p — > 3d tran- 
sitions of iron in the range 16.4-17.2 A. The sample, heated by the thermal 
radiation of a gold spherical hohlraum, was irradiated by the laser ASTERIX 
IV. The plasma is assumed to be in local thermodynamic equilibrium (LTE) 
at a temperature T = 20 eV and a density p = 0.004 g/cm^. An interpre- 
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tation using a Detailed Configuration Accounting (DCA) calculation based 
on SCO code [13] is presented in Fig. [lOl It can be seen that the depar- 
ture from the Gaussian allows one to better reproduce the depth of the 
successive shoulders in the spectrum which correspond to the transitions 
2pi/2 3^3/2, 2p3/2 3(^3/2 and 2^3/2 of several ions. 

The Z-pinch at Sandia National Laboratory was used by Bailey et al. 
[14j to measure absorption of NaBr samples. The main purpose of the 
experiment was to study the 2p — > 3d transitions in bromine ionized into the 
M-shell. Electron temperature and density, obtained from the analysis of the 
sodium lines (which are separated from the bromine lines), are respectively 
50(±4) eV and 3(±1)10^^ cm~^. The spectral resolution is about 1.5 eV. The 
spin-orbit interactions clearly separate the 2pi/2 ~^ 3^3/2 and 2^3/2 3c?5/2 
structures. Fig. [11] displays the experimental spectrum and the calculated 
DCA spectra at T = 47 eV obtained with a Gaussian (dashed curve) or a 
GGu=i function (full line) for the shape of the transition arrays. Some details 
of the 2p 3d structures, hidden in the Gaussian description, appear clearly 
in the new modeling and give a much better agreement with the experiment. 
This provides not only a better identification of the experimental features, 
but also a possible refinement of the temperature and density diagnostic. 

4 The evaluation of high-order moments: a chal- 
lenging task 

4.1 Complexity of the calculation 

Moments /x„ defined in Eq. ^ can be expressed in terms of sums involving 
products of radial integrals. A first difficulty is that the number of terms 
in that sum increases very rapidly as the order increases. For instance, 
assuming that the transition array of interest is characterized by q different 
Slater integrals (R) and r different spin-orbit integrals ((), the number of 
terms of the form R - — R C • • • C is 

n terms p terms 

Nm^^ = Sg{n) X Sr{p), (15) 

where 

fv + t-i\ jv + t-iy. 

^*(") = V J= (t-l)!.! • ^''^ 

Table [3] gives A'^max for each kind of product of such integrals and for 
each of the first four moments n=l, 2, 3 and 4. Obviously, this is a 
maximum since some of the terms do not exist due to angular-momentum 
symmetries. In particular, terms of the kind R - ■ ■ R( containing only one C 
integral are not included since their contribution is zero due to the Lande 
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center-of-gravity rule [U [TT] . One can see that the number of terms increases 
very fast with respect to the order of the moment. The second difficulty 
is that the full derivation of each term of the sum is very cumbersome, 
due to the complex angular-momentum algebra. Uylings has shown [15] 
that the trace of a A;— electron operator in the space of the states of an 

electronic configuration is proportional to ^ ^) ' ^^^'^S second- 

quantization techniques of Judd [4j , one finds that the fourth-order moment 
of — > l^l' involves 1- to 10-electron operators and reads therefore 

-'"n=i;(^;^ (17) 

fc=i ^ ^ 

where coefficients pk have to be determined from special cases and sym- 
metry properties (such as complementarity or anti-complementarity). Par- 
ticular cases are difficult to calculate when dealing with more than two 
electrons (coefficients of fractional parentage are then required to account 
for the Pauli exclusion principle). But the calculation is a hard task even 
with only two electrons ; for instance, considering the simple transition array 
P — > W (corresponding to A^=l), one has: 



ki,k2,k3,k4 

with 



(19) 

where [L] = 2L + 1 and F*=»(//) is a direct Slater integral describing 
electrostatic interactions between two electrons inside orbital /. The quantity 
Aki,k2M,k4 contains also the L-independent product 

[< /||C('=^)||/ >< /||C('=2)||; /||C7(^3)||/ l\\C^''^)\\l >]2, (20) 
where 

<Z||CW||/>=(-l)'[Z](^ J (21) 

is the reduced matrix element of the spherical function operator C^''\ 
with the phase convention of Ref. [TT] . Using the graphical representation 
of 6-j symbols [_16J, a compact graph is obtained for the product of four 
6-j symbols involved in ^A;i,fc2,A:3,fc4 (see Fig. [T2]) . The summation over L 
leads to a 12-j symbol of the second kind, and the summation over L with 
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the phase factor (—1)'^ to a 12- j symbol of the first kind |17j . Therefore, 
the "direct" calculation of the moments is complicated and the final result 
would take many pages. However, some algorithms have been proposed in 
order to evaluate these high-order moments. Karazija et al. [18\ [T9| [20] ex- 
pressed the spectral moments by averages of the products of operators and 
formulated a general group-diagrammatic method for the evaluation of their 
explicit expressions. Oreg et al. [H] used the fact that the moments reduce 
to configuration averages of n-boby symmetrical operators (nBSTOs). For 
that purpose, they introduced the concept of an n-electron minimal con- 
figuration, relative to the actual (N-electron) configuration average. Their 
algorithm uses graphical technique (routine NJGRAF [22j) in order to derive 
the dependence of the averages on the orbital quantum numbers in terms of 
closed diagrams. 



4.2 Estimation of the kurtosis coefficient 

Let us consider the following ideal distribution of line energies and ampli- 
tudes 



Die, a) 



L 



V2 



: exp 



TTV 



e 

2^ 



^exp[-A|a|], 



(22) 



where v is the unweighted variance of the line energies: 

2 



V 



a,b 



a, 6 



ab 



(23) 



e is the line energy, a the line amplitude and L the number of lines. The 
average value of any quantity q{e,a) is given by 



< q >= y J de da D{e,a) q{e,a). 



(24) 



Following Ref. [23], the correlation law between the line energies and 
amplitudes is taken to be 



A2 



7exp[-/3|e|], 



(25) 



where Va is the variance of the line amplitudes, corresponding to the aver- 
age line strength between two levels. The parameters 7 and /? are determined 
by requiring the conservation of the total strength S = J2a b '^ab = L < > 
and of the weighted variance Vw of the line energies: 



S 



1 

52 



a,b 



ab 



L 
S 



(26) 
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The kurtosis is given by 




(27) 



a4 = ^[-2 + (5 + X2)u;], 



(28) 



where uj = v^/v, and X is the root of the fohowing equation: 



(1 + X'^ -a;)exp 



\X^^ 



2 



erfc 



V2 



X 




(29) 



The energy-amphtude correlation is such that, considering the case where 
the spin-orbit interaction is weak, the stronger hues are found cfoser to the 
center of gravity than the weaker ones. Due to these correlations (propen- 
sity rule), the transition array is expected to be sharp. In other words, the 
variance of the line energies is always smaller when it is calculated with a 
weight equal to the line strength (see Eq. (I26p ) than when it is not (see 
Eq. (I23p ). Eq. ()22|) implies that the energy-dependent distribution of line 
strengths S{e) (i.e. the shape of the transition array) can be written 



which means that, if Eq. (j22p is fulfilled, the resulting shape has nec- 
essarily a kurtosis between 3 (Gaussian) and 6 (Laplace). Table H] shows 
the values of estimated from Eq. (f28l) compared to the exact ones (i.e. 
calculated from Cowan's code). 

5 Conclusion 

The UTA formalism gives the first two moments of the distribution of 
electric-dipolar lines. It does not contain any assumption as concerns the 
modeling function. Usually, the chosen distribution is the Gaussian (kur- 
tosis a4=3), but detailed calculations using Cowan's code show that it is 
often not the most proper distribution, since it does not account for the 
skewness and kurtosis coefficients. In this article, several other distribu- 
tions have been studied in an attempt to include the effects of high-order 
moments. Gram-Charlier expansion series, which consists in a Gaussian 
multiplied by a linear combination of Hermite polynomials, can have neg- 
ative values. The Normal Inverse Gaussian function is well suited in order 
to account for asymmetry and flatness. However, the impact of asymmetry 
is usually small in the conditions typical of usual photo-absorption exper- 
iments, and the Generalized Gaussian function seems to be a good choice 




(30) 
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in order to account for the specific effect of the kurtosis [25j. With such a 
modehng, more detailed structures appear in spectra, leading to a better 
agreement with recent 2p — > 3d absorption experiments and improving in 
that way the temperature and density diagnostics. The analytical calcula- 
tion of high-order (n > 3) moments is possible (see the works of Karazija et 
al. [IHIIIHIISO] and Oreg et al. [21J), but difficult. An approximate method, 
relying on a physically realistic distribution of line amplitudes and energies, 
has been proposed in order to estimate the kurtosis. The next step will be 
the calculation of high-order STA (Super Transition Arrays [26j ) moments. 
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Parameter 


Expression 


Role 


a 


3^/3Q4-4a| 
CT(3a4— 5a|^) 


Tail heaviness 




3^3 
(t(3q:4— 5a^) 


Asymmetry 


b 


3(T-y/ 354—503 
334— 4ai^ 


Scale parameter 




Pi '^"•■^ 

-y/3a4— 5a| 


Location 



Table 1: Parameters of the NIG distribution. The quantity a./^ = a4-3 
represents the excess kurtosis. 



order 


exact 


GGjy=1.356 


GC4 


NIG 


Gaussian 


4 


4.15 


4.15 


4.15 


4.15 


3 


6 


33.59 


34.73 


32.22 


38.81 


15 


8 


433.98 


466.46 


346.10 


665.06 


105 


10 


7772 


8954 


4562 


18480 


945 


12 


1.73 10^ 


2.29 10^ 


0.70 10^ 


7.66 10*^ 


10395 


14 


4.52 10** 


7.45 10^ 


1.22 10^ 


4.46 lO*" 


135140 



Table 2: Values of the even reduced centered moments of several distribu- 
tions representing the transition array 3d^4p in Br XII. Exact: calcu- 
lation with Cowan's code; GG: generalized Gaussian function with iv=1.356; 
GC4: fourth-order Gram-Charlier series with 03=0 and a4=4.15; NIG: Nor- 
mal Inverse Gaussian with a=0.82, /3=0, 5=3.17, ;U=154.18. 



Order 


1 


2 


3 


4 


Products involved 


R 


RR, CC 


RRR, RCC, CCC 


RRRR, RRCC, RCCC, CCCC 


Max. number of terms 


4 


16 


54 


150 



Table 3: Number of terms involving products of radial integrals in the first 
moments /i„ of transition array 



Transition array 


kurtosis 


exact 


formula 


Fe V 3^"^ ^ 3£4p 


4.1 


3.8 


Pd VII Ad-^ M-^bp 


4.8 


4.5 


Sn XXIX 3d^ 3dMp 


5.9 


5.1 



Table 4: Estimation of the kurtosis from Eq. (j28p . 
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Figure 1: (Color online) Fourth-order symmetrical Gram-Charlier GC4 dis- 
tributions for several values of the coefficients and 04. The full line 
corresponds to the Gaussian function. 




Figure 2: (Color online) Impact of a skewed Gaussian for modeling the 
opacity lines of a bromine plasma (T=47 eV and /9=0.04 g/cm^) in the 
2p — > 3d range. 
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Figure 3: (Color online) Examples of NIG distributions corresponding to 
different values of the parameters a, /?, 6 and fj,. The full line corresponds 
to tl ^ ■ " 
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Figure 4: (Color online) Validity domain of the NIG distribution (defined 
by jaal < -v/l(a4 — 3)) and positivity domain of the GC4 distribution. 
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Figure 5: (Color online) Examples of GG distributions (corresponding to 
different values of the exponent v). 
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Figure 6: (Color online) Comparison of a transmission spectrum (same as 
Fig. [2]) with a GGu=i profile and a symmetrical GC4 profile both char- 
acterized by the same first four moments (fixing in particular 03=0 and 
04=6). 
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V II 3d^4p-3d^4d 

Foi-ce=6800.04, H|= 4.56 eV, a'=0.46 eV^ a,=1.12, a=9.9&. A=0.2 eV 
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Figure 7: (Color online) Comparison of the line distribution of transition 
array Sd^Ap —>■ Sd^Ad for V II calculated with Cowan's code, and modeled by 
a Gaussian, GG, GC4 and NIG profile. All the distributions are convolved 
by a Gaussian having a FWHM=0.2 eV. 
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Co V 3d'^4s-3d'^4p 

Force=1812.38, (ij = 9.15 eV, 0^=0.52 eV" a3=-0.16, a_,=26.73, A=0.1 eV 
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Figure 8: (Color online) Same as Fig. [7] for the transition array Sd^As 
3d^4p of Co V with a Gaussian linewidth of FWHM=0.1 eV. 
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Figure 9: (Color online) Same as Fig. [7] for the transition array 
3d^4d of Co VI with a Gaussian linewidth of FWHM=0.3 eV. 
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Figure 10: (Color online) Absorption spectrum of iron measured by Chenais- 
Popovics et al. pL2J. The DCA calculations are performed at T=20 eV and 
/?=0.004 g/cm^ assuming either a Gaussian (dashed line) or a GGi^=i (full 
line) profile for the statistical UTA broadening. 
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Figure 11: (Color online) Same as Fig. [10] for the NaBr (in the Br range) 
experiment of Bailey et al. [14J at T=47 eV and p=0.04 g/cm^. 
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Figure 12: Complexity of the calculation of a simple part of the fourth 
moment 114 of the simple transition array P — > II' illustrated by graphical 
methods. 
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